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INTBODUCnON 


Ih«  Ixifrared  spectra  of  crystalline  CS4  and  CD4  In  all  phases  have  recently 
been  studied  by  Savitsky  and  Bomig^.  In  the  low-ten^zature  phase  II  spectra  of 

1.  Q.  B.  Savitsky  and  D.  F.  Homig,  J.  Chen.  Phys.  26^4  (1962). 

pure  CEt  or  CD4  they  observed  fine  structiure  in  the  bending  region  V4  (three  peaks 
separated  by  two  intervals  of  4  cin~^  each^  the  peek  at  the  lowest  frequency  being 
extremely  iiesk)  and  essentially  a  single  peedc  in  the  stret(diing  region,  vs.  The 
spectra  of  dilute  solid  solutions  of  both  isotopic  species  in  one  another  consist 
of  single  bands  %nth  half  widths  3  The  relevant  spectra  are  reparoduced  in 

Figure  1. 

These  observations  indicate  that  the  splittings  in  the  spectra  of  the  pure 
crystals  arise  fron  dynamic  Intemolecular  interactions  rather  than  from  static 
interactions  which  would  show  up  also  in  dilute  solutions. 

There  are  many  structures  for  methane  crystals  that  woiild  split  an  F  mode  into 
two  or  three  Infrared  active  con^onents.  Tha  and  space  groups  with  2 
molecules  per  unit  cell,  or  the  space  gro\qp  having  the  same  packing  with  one 
molecule  per  unit  cell  were  found  to  be  the  most  probable  structures  since  they 
represent  the  most  efficient  hydrogen  packing. 

All  present  calculations  apply  to  the  D|^  as  well  as  the  structure.  The 
latter  is  less  likely  since  it  would  allow  the  molecular  A  and  E  modes  to  become 
active  in  the  crystal  and  they  have  not  been  observed. 

It  may  seem  that  in  these  structures  both  molecules  in  the  unit  cell  are 
parallel  (see  Figure  2)  thus  yielding  a  structure  with  one  molecule  per  unit  cell, 
namely  However,  small  distortions  are  sufficient  to  retain  the  model  of  two 

molecules  per  unit  cell.  In  the  structure  the  molecular  F  modes  sure  split 
into  doublets  rather  than  into  triplets.  This  fact,  taken  alone,  does  not  invali¬ 
date  the  structure  since  one  peak  in  the  observed  q^ctra  is  so  weidc  as  to  be 


uncertain. 
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ABSTRACT 

The  spectrum  of  the  low-temperature  phase  of  crystalline  methane  has  been 
calculated  using  a  model  of  atom-atom  interactions  consisting  of  repulsive, 
attractive,  and  electrostatic  parts  V^,  V^,  and  Vj^.  Detailed  calculations  were 
carried  out  for  the  d|^  and  structures  previously  suggested  for  phase  II 
of  crystalline  methane.  Using  potentials  derived  from  other  sources,  it  is 
possible  to  fnfer  a  potential  for  methane  which  yields  partial  agreement  between 
the  calculated  and  observed  spectra.  In  accordance  with  observations  it  is 
found  that  the  calculated  coupling  splitting  in  the  stretching  region  and  all 
site  splittings  are  srijall  enough  to  be  included  in  the  experimental  band  widths. 
In  the  bending  region  of  the  pure  crystal  spectrum  the  calculated  spacings 
between  the  lines  in  the  triplet  agree  with  experiment  but  the  model  fails  to 
account  for  the  intensity  pattern  and  also  for  the  frequency  shifts  of  the 
solid  solutions  in  either  region. 


^On  leave  of  absence:  Weizraann  Institute  of  Science,  Rehovoth,  Israel, 


-2- 


Aecordlng  to  group  theory  eadi  structure  predicts  the  same  nunher  of  Infrared 
active  bands  for  both  F  nodes  in  pure  crystals  as  well  as  solid  solutions  (See 
OJetble  I). 

To  reconciliate  the  group  theoretical  predictions  with  e^rinental  results 
it  is  required  that  the  numerical  results  for  the  splittings  should  be  such  that 
the  cov^ling  splittings  in  the  bending  region  are  larger  than  in  the  stretching 
region  and  also  larger  than  all  site  splittings. 

In  this  study  the  site  and  coupling  splittings  are  calculated  for  the 
and  structures.  Ohe  calculations  are  based  on  a  model  of  atom-atom  interactions 
recently  Introduced  by  Dows**®^*. 

2.  0.  A.  Dows,  J.  Chem.  Phys.  1^42  (i960). 

3.  D.  A.  Dows,  J.  Chem.  Phys.  202  (I96I). 

4.  D.  A.  Dows,  J.  Chem.  Phys.  2636  (1962). 

CffiNERlL  CONSIDERATIONS 

In  this  section  the  calculation  of  absorption  frequencies  of  molecular  crystals 
is  outlined  and  applied  to  the  case  of  methane.  Inteznolecular  force  constants  for 
the  vibrational  modes  in  the  crystal  are  derived  for  an  appropriate  Intermolecular 
potential  and  the  secular  equation  for  the  unit  cell  is  solved. 

Ihe  following  simplifications  are  made: 

1.  No  mixing  of  stretching  and  bending  modes,  e.g.  the  hydrogens  move 
along  or  perpendicular  to  the  C-H  bonds. 

2.  No  libra tlons. 

3*  Ihe  methane  molecules  in  the  crystal  have  essentially  T^  symmetry. 

4.  Ihe  tetragonal  unit  cell  is  based  on  a  face-centered  cubic  symmetry. 
Ihese  assunptlons  enable  us  to  solve  the  problem  in  a  first  approximation  with¬ 
out  computational  difficulties.®  At  this  stage  it  does  not  seem  Justified  to 

3.  All  calculations  were  carried  out  without  e:qplieitly  taking  into  account  carbon 
recoil.  In  a  later  section,  hdwever,  the  influence' of ' carbon  recoil  on. the 
calculated  spectrum  is  discussed  qualitatively. 
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refine  the  model  before  more  Information  1b  available  about  forces  acting  In  mole¬ 
cular  crystals.  We  are  currently  Investigating  the  Influence  on  the  spectrum  of  a 
change  In  the  symmetry  of  the  molecules  and  also  of  the  unit  cell. 

Symmetry  Coordinates 

Symmetry  coordlnates,S^^,for  the  doubly  occupied  unit  cell  are  built  from 
normal  coordlnatesi  Q,  of  Individual  methane  molecules.  With  the  assumption  of  "no 
mixing",  molecular  normal  coordinates  are  simply  proportional  to  molecular  symmetry 
coordinates  with  proportionality  factors  Las  for  the  stretching  mode  and  L44  for 
the  bending  mode. 

Results  for  the  singly  occv^led  unit  cell  can  be  obtained  In  a  trivial  vay 
from  the  structure.  Iherefore  the  formalism  will  not  be  worked  out  separately. 

From  group  theoretical  consideratiohs.  it  follows  that  each  F  mode  of  the  Indi¬ 
vidual  molecules  gives  rise  to  three  Infrared  active  modes  in  the  crystal,  as  sum¬ 
marized  In  Ihble  I. 


Table  I 

Correlation  of  the  Isolated  Tetrahedron  and  the  Tetrahedron  in  the  Crystal* 
Molecular  Symmetry  Site  Symmetry  Space  Group  Symmetry 

^d  ^2  ^d 


_ > 

As  (z) 

^  (z)  _ ^ 

Bg  (z) 

Fs 

Bg  (y) - 

(x,y) 

§2  (x)  — 

^Underlined  species  are  infrared  active. 


If  we  denote  the  two  molecules  In  the  unit  cell  as  p  and  q  (see  Figure  2),  then 
the  symmetry  coordinates,  S°^,  for  the  unit  cell  may  formally  be  written  as 


S"(A,)  . 

8”«.)  ■  {Op<'x,y>  -  V'x,y>}^‘‘^* 


(1) 


niere  are  two  sets  of  coordinates:  one  built  froo  the  bending  modes  F,  T  ,  F 

X  ^  z 


and  another  set  from  the  stretching  modes  F„y  F, . 

X  y  z 


In  the  free  molecule  the  F.. 


F  ,  and  F  modes  are  degenerate  but  In  a  non-cubic  crystal  they  may  be  different, 
y  z 

As  Indicated  in  Biuatlons  (l)  the  S^(E)  modes  can  be  built  either  from  Q(F^)  oz* 

from  Q(F  ). 
y 

Shifts  and  Splittings 

In  order  to  tfilculate  the  frequencies  of  the  unit  cell  modes  the  force  constant 
matrix  for  the  unit  cell  Is  Introduced 


Xo  +  f, 


PP 


pq 


I 


V  *  'pp_ 

The  G  matrix  In  this  notation  Is  the  identity  matrix.  Solving  the  secular  equation 
for  X  we  obtain  for  the  F  stretching  and  bending  modes 
KBs)  -  Xo  .  f„tQ(P.)l  *  tp,IS(P,)J 


x(As)  -  Xo  ♦  'ppWir.))  • 

X(E,)  .  Xo  .  f„14(P„Jl  *  B™l«(»x.y» 


(S) 


kb.)  ■  Xo  .  -  V‘'’‘’’x,y'’  • 


Here  Xo  >  for  the  Isolated  molecule  and  thus  contains  all  Intramolecular 

interactions. 


PP 


<5) 


and 


"pq“^ 


(4) 
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V  is  the  Inteinolecular  potential,  which  corresponds  to  the  term  in  the 

general  expression  for  the  potential  energy  given  by  Homlg*.  Ihis  part  of  the 


6.  W.  Vedder  and  D.  F.  Homlg,  Advances  in  Spectroscopy  2,  169  (196I). 


potential  contains  cross  terns  between  internal  coordinates  of  all  molecules. 

The  expressions  (2)  for  \  show,  for  each  crystal  mode,  that  X  depends  on  three 
terms.  Ihese  are  Xo>  the  gas  value  for  that  mode,  the  term  f^^  which  produces  a 
shift  and  finally  a  tenn  f^  which  produces  a  splitting, 

^^coupling  “  ^pq  * 

The  coefficient  f^^  in  Equation  (2)  contains  two  terms.  The  first  (the  second 
derivative  with  respect  to  the  normal  coordinate  Q  of  molecule  p)  is  a  static  per> 
turbatlon,  while  the  second  or  "exchange"  term  (the  derivative  with  respect  to 
coordinate  Q  of  molecule  p  and  of  another  molecule  p'  which  is  translationally 
equivalent  to  p)  is  a  dynamic  perturbation.  In  dilute  solid  solutions  of  one  isoto¬ 
pic  species  in  another  f  „  =  0  and  also  the  second  term  of  f  vanishes,  since  at 

pq  PP 

low  concentrations  all  neighbors  are  of  the  other  species  and  do  not  cause  dynamic 

perturbations.  Hence  the  values  of  X  in  a  solid  solution,  X  ,  are  shifted  with 

88 

respect  to  the  gas  value  by  amounts 


x„-Xo  >  s'v/aaj  • 

It  should  be  noted  that  values  for  dX 
around  the  (non  observable)  value  X 
(6)  as 


(6) 


do  not  occur  around  X„  but  rather 
coupling  ss 


crystal? 


which  is  defined  by  Equations  (3)  and 


For  the  stzvcture  with  one  molecule  in  a  unit  cell  on  the  Dg^  site  we 
expect  two  lines  in  the  solid  solution,  due  to  different  static  shifts  in  the  z  and 
x,y  directions 


Xgg(z)  -  Xo  +  a®V/?Q®(Fj 
Xag(x,y)  -  Xo  +  a*V/aQ*(F^^y)  . 
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In  the  pure  crystal  there  will  now  be  two  lines  due  to  different  static  plus-  dynamic 
shifts  In  the  z  and  x,y  directions 

X(B8)  -  Xo  +  6®V/aQ«{F  )  +  a®V/SQ(F^)aQ'(F  ) 

x(B)  =  Xo  +  a*v/aQ®(F  )  a®v/aQ(F^  > 

where  Q'  indicates  a  normal  coordinate  In  a  neighboring  molecule. 

The  pure  crystal  lines  are  the  lines  X(Ba)  and  X(E^)  In  the  D|^  structure. 

The  solid  solution  spectra  are  Identical  for  both  structures. 

ZirrERMQLECULAB  POTEnTIAL 


In  methane  crystals  the  Intermolecxilar  potential  which  Influences  the  vibra¬ 
tional  frequencies  Is  expected  to  be  mainly  due  to  interactions  between  non-bonded 
hydrogens.  Furthezmore,  Intermolecular  Interactions  between  hydrogens  and  carbons 
may  contribute  considerably  to  all  static  shifts. 

For  this  case  the  Intermolecular  potential  can  be  written  In  the  form 


*  ■  f  "'ll’  *  f  I  ■ 


(6) 


where  R,  the  distance  between  hydrogen  1  of  a  given  molecule  and  hydrogen  J  or 
carbon  k  of  a  neighboring  molecule  are  e]q)za88ed  In  Angstrom  units  and  V  In  ergs. 

The  distances  R  are  determined  by  the  geometry  of  the  crystal  irtilch  Is  schema¬ 
tically  drawn  In  Figure  2.  As  an  example  consider  hydrogen  1  of  molecule  p,  or  pi. 

It  has  four  nearest  hydrogens,  labeled  J  >  1,2, ^,4,  at  distances  R^^  ^  "  2.77A; 

these  are  q4,  q^,  (q-a)3>  (q-b)4.  Furthermore,  there  are  six  next-nearest  hydrogens, 
labeled  J  ■  5^6,7,8,9,10,  at  distances  R. .  R™  ■  3«13A.  Two  of  these  hydrogens 
[q2  and  (q-c)2]  belong  to  type  q  molecules,  while  the  other  four  hydrogens  [(p-i«)2, 
(p-»a)4,  (p-fb)2,  (ptb)3]  belong  to  translatlonally  equivalent  molecules  p'  of  type  p"^. 


7.  This  distinction  which  exists,  of  course,  only  in  the  D|.  or  structures  is 
carried  through  all  the  calculations.  For  the  stru^ure  comspondlng  terms 
Involving  p  and  q  molecules  are  combined. 


The  classification  of  neighbor  hydrogens  is  continued  till  all  hydrogens  (69 
In  total)  within  a  radius  of  6a  from  atom  pi  are  labeled  as  is  shown  in  Figure  2  and 
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Appendlx  III.  Similarly  all  carbon  contacts  within  a  radius  of  6k  from  hydrogen 
pi  (l6  In  total)  were  considered  (see  Appendix 4v).  nie  cut-off  at  6a  for  atom- 
atom  Interactions  contributing  to  the  potential  (8)  Is  somewhat  arbitrary.  How- 
evez*^  during  the  course  of  the  calculations  It  was  seen  that  although  at  separations 
of  6a  the  Interactions  In  some  cases  may  remain  Important  (see  (Dables  II  and  III) 
their  contribution  to  Intermolecular  force  constants  become  quite  small.  In  the 
present  study  Index  1  In  Biuatlon  (8)  goes  from  1  to  Index  J  from  1  to  69  and 
Index  k  from  1  to  l6. 

According  to  Biuatlons  (J)  and  (U)  derivatives  of  the  type  have  to  be 

calculated.  For  hydrogen-hydrogen  Interactions,  for  example,  V  Is  defined  in  Equa- 


(9a) 


(9b) 


tlon  (8)  as  an  explicit  function  of  Hence  we  write  for  any  mode  Q 

jay  pjay  jy  Spj  ?pj 

and  elmllarly  for  . 

Here  and  ^  are  normal  displacements  of  hydrogens  1  or  j.  In  the  approxi¬ 
mation  of  "no  mixing"  Is  along  the  C-H  bonds  for  the  stretching  mode  and  for  the 
bending  mode  p'^ is  perpendicular  to  the  C-H  bond  In  the  proper  plane.  Ihe  factors 
(^p/dS  are  elements  of  the  matrix  [U*]  which  transforms  molecular  symmetry  coordi¬ 
nates,  S,  Into  internal  coordinates  Ar  and  dad  (see  Appendix  I).  Ar  Is  the  change 
of  hhe  C-H  bond  length  In  the  stretching  mode,  d  =  1.093A  Is  the  equilibrium  C-H 
bond  length  and  a  Is  the  tetrahedral  angle.  In  the  calculations  It  Is  convenient 
to  express  Internal  coordinates  in  terms  of  9  =  a/2.  Ihe  factors  'S/dQ  are  elements 
of  bhe  matrix  [L]  idrlch  transforms  molecular  normal  coordinates,  Q.,  Into  symmetry 
cooz*dlnates  S. 

As  mentioned  before  hydrogen-carbon  Interactions  contribute  to  the  shifts. 


This  can  be  expressed  by  terms  similar  to  those  In  Q^uatlon  (9Bt)  where  all  J  are 
replaced  by  k  and  where  V  now  stands  for  the  hydrogen -carbon  potential. 
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The  intexnolecular  potential  in  solid  methane  can  be  separated  Into  three 
contributions.  First  there  Is  a  repulsive  potential  of  the  fozia  •  a  exp  (-bR). 
The  paxaoeters  a  and  b  are  not  known  for  atoo-atom  interactions  In  crystalline 
methane.  may  be  taken  froa  the  function  calculated  by  de  Boer®  for  the 

8.  J.  de  Boer,  Biyslca  £,  363  (1942). 

exchange  repulsion  between  non>bonded  H  atoms  In  hydrogen.  Since  It  Is  derived  for 
a  specific  case  the  de  Boer  function  cannot  be  expected  to  apply  rigorously  In  solid 
methane.  Nevertheless  It  has  been  used  with  partial  success  for  calculating  the 
vibrational  spectra  of  CHaCl®*®,  0^4*,  and  paraffin  crystals.® 

9.  R.  Snyder,  J.  Mol.  Spectroscopy  7,  II6  (1961)  . 

Another  expression  used  in  the  calculations  for  given  by  Muller^®,  is  based  on 

10.  A.  Muller,  Proc.  Roy.  Soc.  (London)  A  I78  ,  227  (1941). 


the  measured  compressibilities  of  solid  paraffins. 


The  repulsive  potential  between  hydrogens  and  carbons  was  estimated  to  be 
nR  e 

•  Vp  )  where  the  C-C  repulsive  potential  was  taken  to  be  equal  to  that 


between  two  Ne  atoms^^  ■  2570  x  lO”^®  exp  (-4,26  R). 


11.  J.  0.  Hlrschfelder,  C.  F.  Curtiss,  and  R.  B.  Bird,  Molecular  !nieory  of  Oases 
and  Liquids  (John  Wiley  and  Sons,  Inc.,  New  York,  1954). 


Using  the  de  Boer  function  one  obtains  for  the  repulsive  potential 
V 

R  R  R 

=  [120  exp  (-3.54  R^j)  +  555  exp  (-3.90  R^j^)]  x  lo'^®  erg  (lOa) 


and  with  the  Muller  function  one  obtains 

V_'  -  [4400  exp(-5.0  R.,)  +  3360  exp(-4.63  R.^)]  x  lo"^®  erg  . 


'R  - - -JJ/  •  - - 

Secondly  there  Is  an  attractive  van  der  Waals  type  Interaction  of  the  form 

V  -  +  V^®  ■  -c  R"*  -  c*R“* 

''a  \  ®  "ij  °  ”ik  * 


(10b) 


'9' 


The  coefficients  c  and  c'  were  taken  from  a  calculation  by  Pltcer  and  Catalano^^ 

12.  K.  S.  Fitter  and  E.  Catalano,  J.  Am.  Chen.  Soc.  76,  4844  (1996). 

based  on  the  dispersion  force  formulation  of  constituent  atoms  In  hydrocarbons 

-  [-3.42  r“J  -  8.60  R*J)  X  lo"^®  (erg).  (11) 

A  value  of  c,  about  twice  as  large,  is  found  from  calculation  of  the  dispersion 
Interaction  between  two  H  atoms^^  but  It  should  be  noted  that  this  Is  quite  a 
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are  about  0.12  e  and  about  0.06  e  for  vs  and  V4  reepectlvely^  vbere  e  Is  the  elec¬ 
tronic  charge.  It  is  often  observed  that  in  going  from  the  gas  phase  to  the  solid 
phase  the  intensities  of  infrared  absorptions  Increase,  particularly  for  the 
bending  modes.  Ibis  has  been  demonstrated  for  the  case  of  bentene^^  and  ethylene^^. 


18.  J.  L.  Hollenberg  and  D.  A.  Dowa,  J.  Chem.  Pbys.,  In  press. 

19.  0.  M.  Wleder  and  D.  A.  Dows,  J.  Chem.  Phys.  2990  (1962). 


While  Savitsky  and  Bomlg^  did  not  measure  absolute  Intensities,  F,  of  the  Infrared 
bands  In  solid  methane  it  can  be  seen  from  their  spectra  that  the  ratio  of  Inten¬ 
sities,  r(v4)/r(v3),  measured  on  the  same  samples  is  larger  than  the  ratio  for 
Intensities  in  the  gas  phase.  Ihe  value  e^^^  >  0.1  e  was  Initially  chosen  for  both 
modes  In  methane.  Hence 


E  ''e 


-  [(0.1  e)*R'J  -  (0.1  e)(0.U  e)RjJ]10*»  (erg)^  (12) 

In  addition  calculations  were  made  with  e^^j,  *  0.07  e  and  also  with  0.l4  e  to 

compensate  for  the  lack  of  an  experimental  value  for  e^^^  in  the  solid  state. 

Die  first  and  second  derivatives  with  respect  to  R  of  the  specific  potentials 
(10,  11,  12)  are  numerical  constants  depending  only  on  R.  They  are  listed  In 
Tables  II  and  III  for  all  interatoiolc  contacts  In  crystalline  methane  for  which 


R  <  6a. 

CALCULATIONS  AND  RESULTS 


Ihe  numerical  values  for  the  Individual  factors  In  Qiuations  (9)  will  now  be 
discussed  In  detail, 
a)  The  factors  9p/dS. 

We  choose  a  coordinate  system  for  the  methane  molecule  such  that  the  x,  y,  and  z 
...coordinates  coincide  with  the  three  2-fold  axes  in  the  molecule.  In  the  approx¬ 
imation  of  "no  mixing"  the  molecular  symmetry  coordinates  az« 


-11- 


S®^(A)  ■  (Ari  +  Ara  +  Ara  +  Ar4)2“^ 

S^^(A)  =  d(toia  +  dOja  +  Z^j,4  +  toaa  +  ^034  +  ^034)  6"^^® 

S  (Eg^)  »  d(2daia  -  AQia  -  A0!i4  -  i^aa  -  taa*  +  2^034)12"^'® 

S^*(Ejj)  =  d(-AOi3  +  Aai4  +  Attaa  -  ^024)2"^ 

S®*(F^)  »  (Ari  -  dra  +  dra  -  dr4)2“^ 

S®*(F  )  =  (dri  -  dra  -  dra  +  Lt4)2~^ 

«/ 

S°*(F^)  »  (dri  +  dra  -  dra  -  dr4)2"^ 

S^*^(F^)  =  d(dai3  -  daa4)2"^/® 

s^^(Fy)  =  d(dai4  -  da23)2“^/® 

S^^(F^)  =  d(daia  -  da34)2‘i/® 

The  fhctore  ap/as  are  elements  of  the  matrix  [U*]  rtilch  is  given  in  Appendix  I. 
b)  Ihe  factors  3S/3Q. 

Ihe  factors  as/aQ  are  elements  of  the  matrix  [L].  In  the  approximation  of  no 
mixing 

L  =  (X/F)^/®  .  G^®  ,  (13) 

where  F  and  G  have  the  usual  meaning  in  the  Wilson  notation. ®°  Uhls  results  in 

20.  E.  B.  Wilson,  J.  C.  Declus,  and  P.  C.  Cross,  Molecular  Vibrations,  McGxav- 

Hill,  Inc.,  Nev  York,  1955. 

Las  =  *  1.050  (amu)"^^® 

L44  =  (SUg  +  =  1. 558  (amu)"^^®  . 

More  accurate  values  may  be  obtained  by  inverting  the  complete  matrix  [L**^]  given 
by  Mills®^.  This  results  in 

L.  I.  M.  Mills,  Mol.  Phys.  1,  10?  (1958). 

_ '  <i3 

I 

83  '  1.0497 

^^^2  1 


84  i  -0.1612 


1.5499 
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The  "no  mixing"  assumption,  seems  to  be  Justified  by  the  smallness  of  the  off- 
diagonal  terms  and  by  the  numerical  agreement  between  the  diagonal  elements  in  the 
[L]  matrix  and  the  values  for  Laa  and  L44  obtained  with  Biuation  (13). 
c)  The  derivatives  of  R^j  and 

Calculations  of  the  derivatives  oR/Sp,  5®R/dp®  and  ()®R/d{>j|^?'P j  involved  most 
of  the  computational  work.  The  expressions  for  the  distances  R^^  between  atom  pi 
and  69  neighboring  hydrogens  were  differentiated  with  respect  to 
ASia,  A0I2  (or  ^0^4),  AGla  (or  A0i4)  resulting  in  Ul4  values  for  SR^j/Sp  and 
also  4l4  values  for  the  second  derivatives  which  are  collected  in  Appendices  III 
and  IV  respectively.^® 


22.  Ar.  and  A0.d  are  normal  displacements  of  atom  pi,  while  Ar.  and  A0'.d  are 
normal  displacements  of  hydrogens  in  neighboring  molecules. 


Ihe  carbon  neighbors  contribute  mainly  to  the  static  shifts.  Hence  the  ex¬ 
pressions  for  the  distances  R^j^  between  hydrogen  pi  and  I6  neighboring  carbons 
were  differentiated  with  respect  to  Ar^^,  A0ia,  and  A0i3  resulting  in  48  values 
for  also  48  values  for  These  are  collected  in  Appendix  V. 

A  convenient  way^  to  calculate  derivatives  of  this  type  is  through  the 
following  trigonometric  relations 


.  —  -T _ )  i  ^  i 


(14) 


and 


sin(R''  .  ,p^,  )sin(R,  ,p  .  )cosv 


op  .  . 

'  J 


JO 


(p’^  X  » py 

■  IpJ  IpjI  ’ 

where  y  is  the  dihedral  angle  formed  by  the  three  vectors 


(15) 


\y  Pi. 


-V 

and  pj. 


The  direction  cosines  of  R  can  be  expressed  straightforwardly  in  terms  of 

^  J 

the  parameters  h,  dj  and  0  (see  Figure  2).  The  direction  cosines  of  p  are 


elements  of  the  molecular  matrix  [D]  which  transforms  cartesian  coordinates  into 
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intemal  coordinates  and  which  is  given  in  Appendix  II. 
d)  Results 

We  now  substitute  into  Equations  (9)  the  values  obtained  above  for  the  various 
factors  (except  c>S/?)Q  =  L)  and  take  the  appropriate  sums  over  all  contacts  as 
listed  in  Appendices  VI  -  IX. 

Intermole cular  force  constants  for  H-H  and  H-C  interactions  obtained  with 
the  three  types  of  potentials  considered  in  the  foregoing  section  are  listed  in 
Table  IV.  The  terms  involving  ?iV/3R  and  ?i^V/?R^  are  listed  separately  in  oi-der 
to  show  that  they  are  of  about  equal  Importance.  For  the  repulsive  potential  (10) 
it  may  seem  that  the  terms  with  ('iV/DR  are  smaller  by  a  factor  bR  ~  10.  For  this 
and  other  reasons  terms  with  dV/AR  have  not  been  included  in  previous  calcu¬ 
lations,^"®’'*'^ 

Results  for  the  shifts  and  splittings  (in  cm”^)  are  obtained  from  the  force 
constants  f  through  the  relation 
Av  =  fL^/(8n^c®uo )  • 

In  the  stretching  region  (ua  -  3019  cm""*)  100  dyne  cra”^  is  equivalent  to  O.51 
while  in  the  bending  region  (v4  =  1506  cro”^)  100  dyne  cro"^  is  equivalent  to 
1 . 5O  cm 

From  Equation  (I6)  it  follows  that  for  a  given  value  of  f  the  resulting 
value  of  Au  in  the  bending  region  will  be  larger  than  in  the  stretching  region  by 
a  trivial  factor  (F33G44/F44G33)^^^  which  for  V4  and  03  of  methane  equals  about  5* 

Final  results  for  calculated  shifts  and  splittings  are  listed  in  Table  V. 

Here  we  considered  all  combinations  of  the  three  types  of  interraolecular  poten¬ 
tials  given  above.  For  the  repulsive  potential  these  are;  the  deBoer  potential 
Vj^  and  the  Muller  potential  V^^,.  For  the  electrostatic  potential  these  are  V^, 

2Vj,  and  For  completeness  the  attmetive  potential  was  varied  betveen^^ 

V^^  and  2V^  although  the  latter  value^®  was  derived  for  a  considerably  different 
system. 

The  sum  of  all  attractive  and  repulsive  atom-atom  intemetions  may  be 
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regarded  as  a  measure  of  the  lattice  energy.  This  can  be  compared  to  the  known 
value  of  the  heat  of  sublimation  for  crystalline  methane  which  Is®® 


23.  Handbook  of  Chemistry  and  Physics ^Chemical  Rubber  Publishing  Co., Cleveland,  0. 

2240  cal/mole  =  1700  x  lo"^®  erg/raolecule. 

Good  agreement  is  obtained  for  the  combinations  (V^  +  V^)  and  (V^^,  + 
yielding  respectively  I660  and  20&0  x  lo”^®  erg/molecule.  The  combinations  in¬ 
volving  2V^  yield  values  for  the  lattice  energy  which  are  too  large  by  a  factor  2 
or  more.  Due  to  symmetry  considerations  the  electrostatic  potential  does  not 
contribute  to  the  lattice  energy.  This  was  checked  numerically  for  the  12 
neighbors  arranged  cubically  around  molecule  p  and  indeed  it  was  found  that 

^  . 

E  E  E 

DISCUSSION  OF  RESULTS 

We  now  summarize  the  results  obtained  with  the  present  calculations  and 
compare  them  with  experiment. 

In  good  agreement  with  the  observations  are  the  following  results: 

1)  In  the  bending  region  the  calculated  site  splitting  for  all  cases  (except 
7)  are  small  enough  so  they  may  very  well  be  included  in  the  band  widths. 
In  the  stretching  region  the  site  splitting  is  identically  zero. 

2)  The  coupling  splitting  in  the  stretching  region  in  all  cases  is  also 
small  enou^  to  be  accounted  for  by  the  measured  band  widths. 

3)  The  coupling  splitting  in  the  bending  region  is  of  about  the  right  magni¬ 
tude  for  cases  5 ^  6  or  10  if  we  assume  two  molecules  per  unit  cell.  For 
one  inolecule  per  unit  cell  the  additional  combinations  2  and  12  become 
applicable . 

However,  the  positions  of  the  two  components  Ba  and  E^  are  not  in 
agreement  with  the  previous  assignment^.  Ihe  latter  was  based  on  in¬ 
tensity  considerations  in  the  observed  spectrum,  according  to  which  Bg 
should  be  on  the  low-frequency  side  of  . 
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The  partial  agreement  between  expeilraental  and  calculated  results  seems  to 
verify  the  basic  assumptions  made  at  the  start  of  this  work  about  the  packing  in 
methane  crystals  and  about  tlie  possibility  of  interpreting  the  observed  spectra 
using  a  model  of  atom-atom  interactions. 

The  present  treatment  does  not  allow  us  to  draw  definite  conclusions  about 
the  parameters  in  the  interraolecular  potential.  Out  of  the  combinations  2,  6, 

10,  and  12  which  yield  reasonable  splittings  for  CH4,  5  gives  the  best  result 
for  the  lattice  energy.  The  fact  that  the  calculated  site  splitting  for  this  case 
is  only  barely  acceptable  can  easily  be  rectified  by  increasing  somewhat. 

Very  little  can  be  said  at  this  stage  about  the  shifts.  Both  the  observed 
and  calculated  solid  solution  shifts  for  CH4  show  (for  cases  1,  3^  1>  9> 

11)  that  the  shifts  in  the  stretching  region  are  about  6  cm"^  toward  the  red 
compared  to  the  corresponding  shifts  in  the  bending  region  (see  Table  V).  The 
experimental  shifts  may  be  obtained  by  adding  to  the  calculated  values,  red 
shifts  of  about  8  cm"^  for  the  stretching  as  well  as  the  bending  mode.  These  may 
be  due  to  various  factors  such  as  carbon-carbon  interactions,  the  dielectric 
medium  or  anharraonic  effects  which  are  poorly  known. 

21.  D.  F.  Ball  and  D.  C.  McKean,  Spectrochlm.  Acta  ]^,1029  (1962). 

Deutero-me  thane 

We  now  apply  the  present  calculations  to  the  case  of  deutero-methane .  In 
comparing  CH4  with  CD4  we  first  assume  the  crystal  structure,  the  force  field 
and  the  geometrical  factors  to  be  identical.  The  different  masses  of  H  and  D 
give  rise  to  different  gas  frequencies,  vof  &nd  also  to  different  elements  in 
the  matrix  [L],  which  in  turn  enter  the  results  through  Equation  (16).  If  we 
neglect  the  carbon  movements  in  the  lattice,  then  the  matrices  [D]  for  the  isotopic 
species  are  identical.  In  this  approximation  all  calculated  shifts  and  splittings 
for  CD4  will  be  smaller  than  for  CH4  by  a  factor 
(l2^/l2j^)(vo^/vo®)  =  • 


(17) 
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For  the  bending  region  this  results  in  a  23  °/o  decrease  of  all  calculated 
values.  For  the  solid  solution  frequency  of  CD4  this  decrease  amounts  to  an 
extra  red  shift  of  about  2  cm"^  relative  to  the  solid  solution  shift  of  CH4 
(Table  V;  cases  1,  3,  5^  7,  9^  11)*  Experimentally  we  observe  for  CD4  a  red 
shift  of  3  cm”^  compared  to  a  blue  shift  of  0.5  cra”^  for  CH4.  Half  of  this 
difference  is  accounted  for.  There  is  another  effect  which  produces  a  relative 
red  shift  of  CD4.  This  is  the  weak  coupling  between,  for  example,  the  fre¬ 
quencies  of  CH4  and  the  frequencies  of  CD4  of  the  host  lattice.  This  will  push 
the  energy  levels  of  CH4  to  higher  energies.  Similarly  the  levels  of  CD4  in  a 
host  lattice  of  CH4  will  be  pushed  to  lower  energies. 

A  simple  model  of  one  C-H  bond  in  between  C-D  bonds  or  vice  versa  results  in: 

r  1  (18) 

where  denotes  the  frequency  without  this  particular  effect  and 
the  coupling  splitting  defined  in  Equation  (5).  Substituting  in  (I8)  the  ob¬ 
served  values  for  ^^coupling  shifts  in  the  bending  region  of  -0.05  cm"*^ 

and  0.16  cm  ^  for  CH4  and  CD4  respectively  which  account  for  only  a  fraction  of 
the  difference  in  the  measured  shifts. 

In  the  stretching  rct^ion  the  difference  between  the  observed  solid  solution 
shifts  of  Cn4  and  CD4  cannot  be  explained  by  Equations  (I7,l8). 

The  observed  splittings  in  the  CD4  bending  i-egion  are  about  15  °/o  larger 
than  in  the  CH4  spectrum.  In  order  to  obtain  this  result  the  intermolecular 
force  constants  associated  with  the  splitting  in  the  CD4  spectrum  must  be  about 
50  °/o  larger  than  the  corresponding  force  constants  in  CH4.  This  dlffei'ence 
must  be  caused  by  mass  effects. 

The  first  mass  effect  we  consider  here  involves  the  movements  of  carbon 
atom3--up  till  now  neglected--which  are  relatively  more  important  in  CD4  than 
in  CH4.  A  consequence  of  the  carbon  recoil  during  a  rionnal  vibration  is  a  de¬ 
crease  of  the  movements  of  hydrogens  with  respect  to  a  coordinate  system  fixed 
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in  the  crystal.  For  the  in-phase  crystal  modes  the  carbon  movements  will  not 
influence  intermolecular  interactions  and  hence  will  not  affect  the  spacing 
between  and  Therefore  this  qualitative  explanation  fails  to  account  for 

the  increased  splitting  observed  in  the  CD4  spectrum.  (For  the  out-of-phase 
crystal  modes  recoil  will  cause  a  decrease  of  the  atom-atom  interactions  and 
consequently  a  decrease  of  tlie  spacing  between  E_^  and  E  ) . 

Another  mass  effect  is  the  change  in  the  crystal  dimensions  because  of  the 
difference  in  the  zero  point  amplitudes  of  CH4  and  CD4.  Also  one  might  consider 
the  change  in  each  C-H  bond  length  upon  isotopic  substitution.  This  arises  from 
the  difference  in  the  zero-point  energy  levels  of  CH4  and  CD4  together  with  the 
anharmonicity  in  the  (common)  intramolecular  potential  curve.  However,  the  latter 
two  mass  effects  are  kno\-m  to  be  small  and  it  does  not  seem  likely  that  they  are 
respo’  i.e  for  the  increase  of  50  °/o  in  the  intermolecular  force  constants  In 
going  from  CH4  to  CD4. 
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TABLE  II 


Values  for  the  First  and  Second  Derivatives  of  the  Hydrogen-Hydrogen 
Potential  in  Methane  Crystals  (Units:  dyne  x  lO”®  and  dyne  cm  ^  Respectively) 
Repulsive  Potential  Attractive  Electrostatic 


de  Boer  Muller  Potential  Potential 


■'ij 

^  , 

^  i 

'■"ij 

^«"lj 

2.77 

-235 

852 

-216 

1079  ■ 

164 

-4i6 

-301 

217 

3.13 

-65 

229 

-35 

173 

69 

-155 

-235 

150 

4.06 

-3 

9 

0 

2 

11 

-20 

-l4o 

69 

4.13 

-2 

7 

0 

1 

10 

-17 

-135 

65 

4.25 

-1 

i 

k 

8 

-13 

-128 

1 

60 

4.50 

1 

2 

! 

6 

-9 

i  -114 

1 

51 

4.74 

'  0 

1 

1 

4 

-0 

1  -103 

44 

5.18 

2 

-3 

1  -86 

33 

5.57 

i 

1 

-2 

1  -75 

27 

5.60 

1 

-1 

1  -69 

24 

» 


TABLE  III 


Values  for  the  First  and  Second  Derivatives  of  the  Hydrogen-Carbon  Potential 
in  Methane  Crystals  (Units:  dyne  x  lo”®  and  dyne  cm"^  Respectively) 


TABLE  IV  A 


Intermolecular  Force  Constants  in  Solid  Methane*  (Units:  dyne  cm”^) 


i 

j  Repulsive  Forces 

Attractive 

Electrostatic 

•  Total 

Bending  Region 

1 

de  Boer 

1 

Muller 

Forces 

Forces 

a)  de  Boer 

1 

1 

1 

1 

b)  Muller 

Solid  solution  shift 

1356+92 

l400+60 

-900-220 

1398-1210 

a)  4l4 

-343-5^ 

-281-39 

302+114 

-1434+1313 

z 

1051 

ii4o 

-704 

07 

0)  505 

Solid  solution  shift 

1613+112 

1725+68  1 

-1029-263 

1467-1255 

a)  560 

-317-  53 

-258-40 

284+108 

-1399+1292 

b)  700 

y 

1355 

14$5 

-900 

105 

Crystal  shift 

1149+92 

1244+60 

-749-220 

1155-1210 

a)  468 

-317-54 

-158-39 

262+114 

-1067+1313 

z 

bio 

1107 

-593 

191 

b)  705 

Crystal  shift 

1116+112 

1357+68 

-655-265 

865-1255 

a)  51 

-310-  53 

-255-40 

253+108 

-1159+1292 

X/  y 

I  005 

“ll^O 

-55Y 

- - 

b)  516 

Coupling  splitting 

1  -425 

-381 

305 

1 

CD 

a)  -5 

z 

206 

-219 

1 

-122 

-191 

-155 

150 

1 

b)  23 

Coupling  splitting 

1  -208 

-275 

156 

-381 

a)  242 

X,  y 

182 

!  151 
'  -122 

1 

1 

-156 

0 

649 

b)  l46 

*See  footnote  in  TABLE  IV  B 


TABLE  IV  B 


Intermolecular  Force  Constants  in  Solid  Methane*  (Units:  dyne  cm”^) 


Repulsive  Forces 

Attractive 

Electrostatic 

Total 

Stretching  Region 

de  Boer 

Muller 

Forces 

Forces 

a)  de  Boer 

1 

1 

1 

b)  Muller 

Solid  solution  shift  ' 

1790+5^8 

2135+502 

,  -1055-502 

1391-1550 

a)  737 

1 

1 

-294-  15 

-227-  8  1 

280+  71 

-1436+1307 

b)  1110 

X,  y,  z  , 

202$ 

240^  i 

!  -1204 

-tJtJ 

Crystal  shift  ! 

1788+548 

2135+502 

-1049-502 

1347-1350 

a)  922 

i  -293-  15 

-227-  8 

268+  71 

-1198+1307 

b)  1296 

z 

2028 

2402 

-1212 

ioT^ 

Crystal  shift 

1  1735+548 

2105+502 

1  -984-502 

1161-1350 

a)  669 

-285-  15 

-222-  8 

1  263+  71 

-1282+1307 

b)  1061 

X,  y 

1985  1 

■^377  : 

■  ^1155 

-164 

Coupling  splitting 

372 

-400 

243 

-576 

a)  -276 

68 

68 

-51 

212 

z 

-3oir 

-352 

b)-304 

Coupling  splitting 

74 

1 

i  67 

-16 

-149 

a)  190 

X,  y 

25 

99 

1  12. 

1  79 

1 

1 

-39 

-55 

295 

b)  170 

! 

1 

1 

*The  four  entries  are  arranged  as  follows 
Upper  left:  (DR^j/dS)^ 

Lower  left: 

Upper  right :(3®V^*^/^R^j^)  (aR^^^/QS)^ 
Lower  right: 


TABLE  V 


Shifts  and  Splittings  for  Various  Combinations  of  the  Three  Intermolecular 
Interactions  and  in  Solid  Methane  (All  Entries  In  cm“^) 


1 

Interaction 

BENDING 

R  E  G  I 

0  N  j 

STRETCHING  REGION 

Solid 

Solution 

Shift 

i 

(a^iy+2:)/2! 

Site 

Splitting 

x,y-z 

Coupling  Splitting 

Solid 

Solution 

Shift 

(x,y+z)fe 

Coupling  Splitting 

Ba-E^ 

+ 

1 

1 

Ba-E^ 

E_^-E_ 

1 

- 1 

V  +v  +v 

V  E 

1 

7.6 

t  2.2 

2.6 

! 

7.5 

HI 

wm 

1.1 

2 

V  +2V  +V 

''r  ^  V  E 

-*+.9 

-0.7 

4.4 

7.5 

B 

0.8 

3 

V  +V  +2V 

R  A  ^E 

8.9 

2.8 

6.5 

15.9 

2.0 

-0.7 

2.0 

k 

V  +2V  +2V 

R  A  ^  E 

-3.5 

1 

i  -0.1 

8.2 

15.9 

-1.7 

-0.2 

1.7 

5 

9.4 

3.0 

4.1 

4.5 

I  3.4 

1 

-0.7 

1.0 

6 

-3.1 

1 

0.0 

1 

5.9 

4.5 

-0.3 

-0.1 

0.7 

7 

10.7 

3.6 

6.4 

12.9 

3.1 

-0.8 

1.9 

8 

Vj,,t2V^.2V^ 

-1.0 

1 

0.6 

9.7 

12.9 

-0.5 

-0.2 

1  1.6 

9 

vw= 

i  6.9 

1.9 

0.8 

3.3 

2.4 

-0.6 

0.7 

10 

1 

i  1 

-5.5 

2.5 

3.3 

-1.3 

0.0 

0.3 

11 

1  8*7 

!  2.3 

0.3 

3.5 

-0.6 

0.6 

12 

■I 

1  -3.8 

4.2 

0.3 

-0.1 

1 

i  ..  J 

-0.1 

0.2 

1 

. ■■■) 

MMMi 

■nnnm 

Observed 

1 

0.5 

<3 

-3.9 

<4 

<4 

:d4 

Observed  ; 

I 

-3.9 

1 

i  ^ 

1 

'  -4.5 

B 

-7.-2 

<4 

i 

<4 

1 

APPENDIX  I 


The  [U']  Matrix  for  Molecular  Methane 

f(Ar)  S(pf )  S(F^^)  S(Ff )  S(F^^)  S(fJ‘^)  S(Ff )  S(Ef )  S(E^^)  SCAi"^) 

An  !  2"^  2’^  2"^  2‘^ 

Are  ;  2“^  -2“^  -2“^  2"^ 

Ara  2“^  2"^  -2“^  -2"^ 

Ar4  2"^  -2"^  2"^  -2“^ 


1 

Aotiad 

1 

0 

0 

2-1/2  5-1/2  0 

6-1/2 

Aa34<i 

0 

0 

.2-1/2  5-1/2  0 

6-1/2 

AOiad  : 

0 

0  -12'^/®  -2"^ 

6-1/2 

Aaa4d  1 

.2-1/2 

0 

0  -12’^/^  -2"^ 

6-1/2 

Aaa4d  1 

0 

2-1/2 

0  -12’’^/®  2"^ 

6-1/2 

1 

Aotaad  i 

0 

_2“i/2 

0  -12“^/®  2"^ 

6-1/2 

APPERDIX  U 


Ihe  [D]  Matrix  for  H61eeular  Matbane''^ 


Internal 

Coordi¬ 

nates 

Nomall- 

sation 

Factor 

dxi 

dyi 

dzx 

dxa 

dya 

dze 

dxs 

dya 

dZs 

dXi 

dye 

dZ4 

dn 

j-i/a 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

dra 

3-1/a 

0 

0 

0 

-1 

-1 

1 

0 

0 

0 

0 

0 

0 

dra 

3-1/a 

0 

0 

0 

0 

0 

0 

1 

.1 

-1 

0 

0 

0 

dr4 

1 

10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

1 

-1 

daiad 

6-i»a 

1 

1 

-2 

-1 

-1 

-2 

0 

0 

0 

0 

0 

0 

dCitS4d 

g-i/a 

0 

0 

0 

0 

0 

0 

1 

-1 

2 

-1 

1 

X 

dcxisd 

g-i/a 

-2 

1 

1 

0 

0 

0 

-2 

-1 

-1 

0 

0 

0 

dOae^ 

5-1/a 

0 

0 

0 

2 

-1 

1 

0 

0 

0 

2 

1 

-1 

da24d 

•* 

5-1/a 

1 

-2 

1 

0 

0 

0 

0 

0 

0 

-1 

-2 

-1 

dOfes^ 

6-1/a 

0 

0 

0 

-1 

2 

1 

1 

2 

-1 

0 

0 

0 

*Slx  rows  for  the  conditions  of  no  translation  and  rotation  aure  omitted;  so  are  the 
columns  for  the  three  coordlxiates  of  the  ccuihon  atom  and  one  for  the  redundancy 
condition. 


APPENDIX  III 


First  Derivatives  of  the  Intermolecular  Hydrogen  Distances  in  Crystalline  Methane* 


Neighbor 

Hydrogen 

J 

'’ij 

ar.  5^ 

3R 

ijj-i 

dR,  . 

aR. . 

SR 

q4 

1 

2.77 

'-0.s^44  -0.083 

0.066 

-0.312 

0.297 

qj 

2 

2.77 

-0.944  -0.083 

0.066 

0.246 

0.675 

-0.972 

(q-a)5 

5 

2.77 

-0.083  -0.944 

0.675 

-0.972 

0.066 

0.246 

(q-b)4 

4 

2.77 

-0.083  -0.944 

0.675 

0.297 

0.066 

-0.312 

q2 

5 

3.13 

-0.834  0.242 

0.552 

-0.276 

-0.970 

•  0.485 

(q-c)2 

6 

3.15 

0.242  -0.834 

-0.970 

0.485 

0.552 

-0.276 

(p+a)2 

7 

3.13 

-0.296  -0.296 

-0.209 

0.912 

-0.209 

0.912 

(p+a)4 

8 

3.13 

-0.296  -0.296 

-0.703 

0.912 

-0.703 

0.912 

(p+b)2 

9 

3.13 

-0.296  -0.296 

-0.209 

-0.705 

-0.209 

-0.703 

(p+b)5 

].0 

3.13 

-0.296  -0.296 

-0.703 

-0.209 

-0.703 

-0.209 

(p+a+b)2 

11 

4.o6 

-0.8i6  -O.O16 

-0.577 

0.289 

-0.577 

0.289 

qi 

12 

4.13 

-0,986  0.986 

0.169 

-0.085 

-0.169 

0.085 

(q-c)l 

13 

4.13 

-0.169  0.169 

-0.986 

0.495 

0.986 

-0.493 

(q-h)l 

14 

4.13 

-0.4ob  0.4o8 

0.577 

0.324 

-0.577 

-0.324 

(q-b-c)l 

15 

^.13 

0.4o8  -0.408 

-0.577 

0.901 

0.577 

-0.901 

(q-a-b)l 

l6 

4.13 

0.169  -0.169 

0.986 

-0,493 

-0.986 

0.493 

(q-a-b-c)l 

1? 

4.13 

0.986  -0.986 

-0.169 

0.085 

0.169 

-0.085 

(q-a)i 

18 

4.13 

-0.4o8  o.4o8 

0.577 

-0.901 

-0.577 

0.901 

(q-a-c)l 

19 

4.13 

0.4o8  -o.4o8 

-0.577 

-0.324 

0.577 

0.324 

(p+a)l 

20 

4.13 

-0.577  0.577 

-0.4o8 

0.816 

0.4o8 

-0.816 

(p-b)l 

21 

4.13 

0.577  -0.577 

0.4o8 

0.4o8 

-0.4o8 

-0.4o8 

(p-a)l 

22 

4.13 

0.577  -0.577 

0.4o8 

-0.816 

-0.4o8 

0.816 

(p+b)l 

23 

4.13 

-0.577  0.577 

-o.4o8 

-o.4o8 

0.4o8 

0.4o8 

APPENDIX  III  (Continued) 

First  Deidvatives  of  the  Interraolecular  Hydrogen  Distances  in  Crystalline  Methane* 


Neighbor 

Hydrogen 

J 

^^i 

’■'ij 

SR.  . 

^0i2 

SR 

--iid“^ 

hQxs 

3H,, 

SR.  . 

^0X4 

(q-b)3 

24 

4.25 

-0.054 

0.507 

0.439 

0.557 

0.836 

-0.237 

(<l-a)4 

25 

4.25 

-0.054 

0.507 

0.439 

-0.996 

0.856 

-0.599 

(q-a-b)3 

26 

4.25 

0.507 

-0.054 

0.836 

-0.237 

0.439 

0.557 

(q-a-b)4 

2? 

4.25 

0.507 

-0.054 

0.836 

-0.599 

0.439 

-0.996 

(q-b)2 

28 

4.50 

-0.051 

0.699 

0.759 

0.182 

-0.301 

0.712 

(q-b-c)2 

29 

4.50 

0.699 

-0.051 

-0.301 

0.712 

0.759 

0.182 

(q-a)2 

30 

4.50 

-0.051 

0.699 

0.759 

-0.941 

-0.301 

-0.411 

(q-a-c)2 

31 

4.50 

0.699 

-0.051 

-0.301 

-0.411 

0.759 

-0.941 

(p+a)5 

32 

4.50 

-0.206 

0.854 

-0.489 

0.979 

0.260 

-0.520 

(p+b)4 

33 

4,50 

-0.206 

0.854 

-0.489 

-0.489 

0.260 

0.260 

(p-a)3 

34 

4.50 

0.854 

-0.206 

0.260 

-0.520 

-0.489 

0.979 

(p-b)4 

35 

4.50 

0.854 

-0.206 

0.260 

0.260 

-0.489 

-0.489 

(q-c)3 

36 

4.73 

0.161 

0.66^^ 

-0.969 

0.647 

-0.613 

-0.065 

(q-c)4 

37 

4.73 

0.161 

0.664 

-0.969 

0.322 

-0.613 

0.678 

(q-b-c)4 

36 

4.73 

0.664 

0.161 

-0.613 

0.678 

-0.969 

0.322 

(q-a-c)3 

39 

4.73 

0.664 

0.161 

-0.613 

-0.065 

-0.969 

0.647 

(p+a+b)3 

4o 

5.16 

-0.639 

0.281 

-0.750 

0.524 

-0.099 

-0.777 

(p+a+b)4 

4l 

5.16 

-0.639 

0.281 

-0.750 

0.225 

-0.099 

0.376 

(p+a-b)4 

42 

5.18 

0.281 

-0.639 

-0.099 

0.876 

-0.750 

0.225 

(p-a+b)3 

43 

5.18 

0.281 

-0.639 

-0.099 

-0.777/ 

-0.750 

0.524 

(p-a)2 

44 

5.54 

0.694 

0.694 

0.491 

-0.702 

0.491 

-0.702 

(p-b)2 

45 

5.54 

0,694 

0.694 

0.491 

0.211 

0.491 

0.211 

(p-a)4 

46 

5.54 

0.694 

0.694 

0.211 

-0.702 

0.211 

-0.702 

(p-b)3 

47 

5.54 

0.694 

0.694 

0.211 

0.491 

0.211 

0.491 

or 


APPENDIX  III  (Continued) 

First  Derivatives  of  the  Intermolecular  Hydrogen  Distances  in  Crystalline  Methane* 


Neighbor 

Hydrogen 

i 

ij  drj^  dr^ 

0012 

SR  _ 

dR,  , 
iJj-i 

dR,  . 

351?* 

1 

P 

1 

O' 

ro 

48 

5.54  0.390  0.998 

0.921 

-0.461 

0.060 

-0.030 

(q-a-b-c)2 

49 

5.54  0.998  0.390 

0.060 

-0.030 

0.921 

-0.461 

(n4a)4 

50 

5.60  -0.893  -0.466 

-0.269 

0.447 

0.033 

0.749 

(q+b)3 

51 

5.60  -0.893  -0.466 

-0.269 

-0.179 

0.033 

-0.782 

(q4a-b)4 

52 

5.60  -0.466  -0.893 

0.033 

0.749 

-0.269 

0.447 

(q-a+b)3 

53 

5.60  -0.466  -0.893 

0.033 

-O.7B2 

-0.269 

-0.179 

(q-b-c)3 

54 

5.73  0.549  0.965 

-0. 506 

0.829 

-0.212 

0.241 

(q-a-b-c)3 

55 

5.73  0.965  0.549 

-0.212 

0.241 

-0.506 

0.829 

(q-n-b-c)4 

56 

5.73  0.965  0.549 

-0.212 

-0.029 

-0.506 

-0.325 

(q-a-c)4 

57 

5.73  0.549  0.965 

-0.506 

-0.323 

-0.212 

-0.029 

(q+b)2 

58 

5.79  -0.863  -0.281 

0.007 

-0.441 

-o.ai6 

-0.029 

(q4a)2 

59 

5.79  -0.863  -0.281 

0.C07 

0.434 

-0.816 

0.645 

{q4b-c)2 

60 

5.79  -0.281  -0.863 

-0.816 

-0.029 

0.007 

-0.441 

(q4a-c)2 

61 

5.79  -0.281  -0.863 

-0.816 

0.845 

0.007 

0.434 

(p4a4b)l 

62 

5.84  -0.816  0.816 

-0. 577 

0.289. 

0.577 

-0.289 

(p4a-b)i 

63 

5.84  0  0 

0 

0.866 

0 

-0.866 

(p-a-b)l 

64 

5.84  0.816  -0.816 

0.577 

-0.289 

-0.577 

0.289 

(p-a4b)l 

65 

5.84  0  0 

0 

-0.866 

0 

0.866 

+(p-c)l 

66 

5.84  0.577  -0.577 

-0.816 

0.4o8 

0.816 

-0.4o8 

4(p4c)l 

61 

5.84  -0.577  0.577 

0.816 

-0.4o8 

-0.816 

0.408 

+(p-tc)5 

68 

4,75  -o.4o3  -0.403 

0.8y6 

-0.285 

0.896 

-0.285 

4(p4c)4 

69 

4.75  -o.4o3  -0.4o3 

0.O96 

-0.611 

0.896 

-0.611 

*Ar^  and  hod  are  noriEal  displacementc  of  atom  pi,  while  Ar  and  AO 'd  are  normal 
displacements  of  hydrogens  in  neighboring  molecules.  ^ 
t-not  indicated  in  Figure  2 


APPENDIX  IV 


Second  Derivatives  of  the  Interraolecular  Hydrogen  Distances  in  Crystalline  Methane* 


Neighbor 

Hydrogen 

J 

dr, df,  dr. 
i  i  j 

defa 

ijj-2 

S013 

i.1  .  j-2  __ 

501 pSOia 

h^E.. 

50135013 

501350  A.-. 

q4 

1 

0.039  0.092 

0.360 

0.326 

0.225 

0.274 

q3 

2 

0.039  0.092 

0.560 

0.539 

0.225 

-0.034 

(d-a)3 

3 

0.339  0.092 

0.197 

0.020 

0.225 

-0.034 

(q-a)4 

4 

0.359  0.092 

0.197 

0.330 

0.225 

0.27^ 

q2 

5 

0.097  0.171 

0.222 

0.295 

0.065 

0.255 

(q-c)2 

6 

0.300  0.171 

0.019 

0.244 

0.065 

0.255 

(p+a)2 

7 

0.291  0.078 

0.305 

0.054 

-0.120 

-0.053 

(lHll)4 

8 

0.291  0.078 

0.162 

0.05^ 

0.055 

-0.053 

(p+b)2 

9 

0.291  0.078 

0.305 

0.162 

-0.120 

0.055 

{p+4)3 

10 

0.291  0,078 

('.162 

0.305 

0.055 

-0.120 

( p+a+b)2 

11 

0.082  -0.082 

0.164 

U.226 

-0. 16^+ 

0.144 

qi 

12 

0.007  -0.007 

0.255 

0.24o 

-0.235 

-0.240 

(q-c)l 

13 

0.235  -0.235 

0.007 

0.183 

-0.007 

-0.183 

(q-b)i 

l4 

0.202  -0.202 

<1.162 

0.217 

-0.162 

-0.217 

(q-b-c)l 

13 

0.202  -0.202 

0.162 

o.o46 

-0.162 

-0.046 

(q-a-b)l 

16 

0.235  -0.235 

0.007 

0.183 

-0.007 

-0.183 

(q-a-b-c)l 

17 

0.007  -0.007 

0.235 

0.24o 

-0.235 

-0.240 

(q-a)i 

lO 

0.202  -0.202 

0.162 

o.o46 

-0.162 

-0.046 

(q-a-c)l 

19 

0.202  -0.202 

0 . 162 

C.217 

-0.162 

-0.21? 

(p+a)l 

?o 

0.162  -0.162 

0,202 

0.081 

-0.202 

-0.081 

(p-b)l 

21 

0.162  -0.162 

0.202 

0.202 

-0.202 

-0.202 

(F-a)l 

O'') 

0.162  -G.I62 

0  •  202 

0.081 

-0.202 

-0.081 

(piV,)l 

23 

0.162  -0.162 

0 . 202 

0.202 

-0.202 

-0.202 

APPENDIX  IV  (Continued) 


Second  Derivatives  of  the  Interraolecuiar  Ilydroscn  Diotarices  in  Crj'^stalline  Methane* 


Neighbor 

Hydrogen 


cOisdSi 


or 


(q-b)3 

24 

0.234 

0.085 

0.190 

(q-a)4 

25 

0.234 

0.085 

0.190 

1 

P 

1 

O' 

26 

0.175 

0.085 

0.071 

(q-a-b)4 

27 

0.175 

0.085 

0.071 

(q-b)2 

28 

0,222 

0.082 

0.094 

(q-b-c)2 

29 

0.il4 

0.082 

0.202 

(q-a)2 

30 

0.222 

0.082 

0.094 

(q-n-c)2 

31 

0.114 

0.082 

0.202 

(p+a)3 

32 

0.213 

0.113 

0.169 

(p+b))| 

33 

0.213 

0.113 

0.169 

(p-a)3 

34 

0.0  Co 

0.115 

0.207 

(p-b)4 

35 

o.oCo 

0,113 

0.207 

(‘l-c)3 

36 

0.206 

0.042 

0.014 

(q-c)4 

37 

0.206 

O.042 

0.014 

(q-b-c)4 

38 

0.118 

0.042 

0.139 

1 

P 

I 

O 

39 

0.116 

0.042 

0.159 

(p+u+b)3 

4o 

0.114 

0.099 

0.084 

(p-ta-hb)4 

4l 

0.1.14 

0.099 

0.084 

(p-m-b)4 

12 

0.178 

0.099 

0.191 

(p-a+b)3 

43 

0.1 70 

0.099 

0.191 

(p-a)2 

14 

0.094 

-0.027 

0.137 

(p-b)2 

45 

0.094 

-0.027 

0.137 

(p-n)4 

46 

0.094 

■0.027 

0.172 

(p-b)3 

47 

0.094 

-0.027 

0.172 

0.162 

O.O  ^0 

-0.o47 

0.002 

0.070 

0.016 

0.222 

0.070 

-0.047 

0.151 

0.070 

0.016 

0.215 

-0.058 

0.119 

0.110 

-0.058 

0.119 

0.026 

-0.058 

0.062 

0. 185 

-0.058 

0.062 

6.009 

0.177 

0.039 

0.169 

0.177 

0.177 

0.162 

0.177 

0.039 

0.207 

0.177 

0.177 

0.123 

0.015 

-0.062 

0.189 

0.015 

0.095 

0.114 

0.015 

0.095 

0.210 

0.015 

-0.062 

0.l40 

0.114 

0.014 

0.IH5 

0.114 

0.091 

0.045 

■  0.114 

0.091 

0.076 

0.1l4 

0.014 

0.092 

-0.103 

0.031 

0.172 

-0.103 

0.112 

i'.092 

0.112 

0.031 

0.1.12 


-0.103 


APPENDIX  IV  (Continued) 

Second  Derivatives  of  the  Inteiroolecular  Hydrogen  Distances  in  Crystalline  Methane* 


Neighbor 

Hydrogen 

J 

1  Jj"2 
6012 

lJj-2 
■■  d 

'^*13 

A 

50X26012 

5^ 

A-S 

601260^4 

6013613 

50 1360^4 

(q-a-b)2 

48 

0.155 

-0.010 

0.027 

0.142 

-0.07c 

0.118 

(q-a-b-c 

)2  49 

0.001 

-0.010 

o.ieo 

0.180 

-0.070 

0.118 

(q+a) 4 

50 

0.056 

-0.015 

0.166 

0.143 

0.121 

0.059 

(q+t)3 

51 

0.036 

-C.015 

0.166 

0.173 

0.121 

0.070 

(q+a-b)4 

52 

o.i4o 

-0.015 

0.178 

0.078 

0.121 

0.059 

(q-a+b)3 

53 

0.140 

-0.015 

0.178 

0.069 

0.121 

0.070 

(q.b-c)3 

54 

0.122 

-0.034 

0.130 

0.054 

0.098 

-0.093 

(q-a-b-c 

)3  55 

0.012 

-0.034 

0.167 

0.073 

0.098 

-0.093 

(q-a-b-c 

)4  56 

0.012 

-0.034 

0.167 

0.174 

0.098 

0.115 

(q-a-c)4 

57 

0.122 

-0.034 

0.130 

0.156 

0.098 

0.115 

( q+b ) 2 

58 

0.044 

0.016 

0.173 

0.159 

-0.056 

0.115 

(q+a)2 

59 

0.044 

0.016 

0.175 

0.140 

-0.056 

0.052 

(q+b-c)2 

6o 

0.159 

0.016 

0.058 

0.173 

-0.056 

0.115 

(q+a-c)2 

6i 

0.159 

0.016 

0.058 

0.049 

-0.056 

0.052 

(p+a+b)l 

62 

0.057 

-0.057 

0.114 

0.157 

-0.114 

-0.157 

(p+a-b)l 

63 

0.171 

-0.171 

0.171 

0.043 

-0.171 

-0.043 

(p-a-b)l 

0.057 

-0.057 

0.114 

0.157 

-0.114 

-0.157 

(p-a+b)l 

65 

0.171 

-0.171 

0.171 

0.043 

-0.171 

-o.o43 

+(p-c)l 

66 

0.114 

-0.114 

0.057 

0.145 

-0.057 

-0.145 

+(p+c)l 

67 

0.114 

-0.114 

0.057 

0.145 

-0.057 

-0.143 

+(p+c)3 

68 

0.176 

0.056 

0.04l 

0.194 

-0.029 

-0.087 

+(p+c)4 

69 

0.176 

0.056 

o.o4i 

0.132 

-0.029 

0.062 

-2 


-2 


*Ar  and  AGd  are  normal  displacements  of  atom  pi,  while  Ar 
displacements  of  hydrogens  in  neiglihorlng  molecules. 

+aoL  indicated  in  Figure  2. 


and  A0’d  are  normal 


or 
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Shifts  and  S|pllttings  for  Various  Conblnatlons  of  the  Diras  IntsriBOlecular 

*  /  -1  V 

Interactions  and  In  Solid  Methane  (All  Ent  ries  In  ea 


Interaction 

— mmw 

Tumo - 

Solid 

Solution 

Shift 

(Xiy+i)/2 

Site 

Splitting 

x,y-2 

Coupling  Splitting 

Solid 

Solution 

Shift 

(x,y+z)/fe 

Coupling  Splitting 

Ba-E^ 

1  1 

E^.E 

T  • 

Ba-E^ 

E^-E_ 

1 

VV^E 

7.6 

2.2 

2.6 

7.5 

2.2 

•0#  6 

1.1 

2 

-4.9 

-0.7 

4.4 

7.5 

-1.4 

0.1 

0,8 

3 

VV^Ve 

8.9 

2.8 

6.5 

15.9 

2.0 

-0.7 

2.0 

4 

-3.5 

-0.1 

8.2 

15.9 

-1.7 

*0»  2 

1.7 

5 

9.4 

3.0 

4.1 

4.5 

3.4 

-0.7 

1.0 

6 

-3.1 

0.0 

5.9 

4.5 

-0.5 

-0,1 

0.7 

7 

10.7 

3.6 

6.4 

12.9 

3.1 

-0.8 

1.9 

8 

-1.8 

0.6 

9.7 

12.9 

-0.5 

-0,2 

1.6 

9 

WV^ 

6.9 

1.9 

0.8 

3.3 

2.4 

-0.6 

0.7 

10 

V»A*V= 

-5.5 

2.5 

3.3 

-1.3 

0.0 

0.3 

11 

8.7 

2.3 

0.3 

3.5 

-0,6 

0.6 

12 

-3.8 

4.2 

0.3 

-0,1 

-0.1 

0,2 

a* 

Observed 

1 

0.5 

<3 

-3.9 

4.1 

-7.6 

<4 

r*"***» 

<4 

■D4 

Observed  | 

I 

-3.9 

i 

<1 

-4.3 

4.9 

-7.2 

<4 

<4 

APPENDIX  VI 


Values  of  for  Crystalline  Methane* 


J 

’*1.1 

Type 

STRETCH 

X,  y,  z 

BEND 

z 

BEND 

X,  y 

1-4 

2.77 

<1 

1.796 

0.920 

1.191 

5-10 

3.13 

p,q 

1.107 

2.322 

2.513 

11 

4.06 

p 

0.666 

0.333 

0.084 

12  -  23 

4.13 

p^q 

4.000 

4.000 

4.000 

24  -  27 

4.25 

q 

0.520 

1.783 

1.717 

ro 

CD 

4.50 

p^q 

2.526 

1.948 

3.129 

36  -  39 

4.73 

q 

0.934 

2.630 

0.987 

1 

o 

5.18 

p 

0.974 

1.146 

1.697 

44  -  49 

5.54 

p>q 

3.066 

1.424 

1.485 

50  -  53 

5.60 

q 

2.028 

0.146 

1.4o5 

54  -  57 

5.73 

q 

2.464 

0.602 

1.573 

58  -  61 

5.79 

q 

1.648 

1.332 

1.097 

62  -  67 

5.84 

p 

2.000 

2.000 

2.000 

68  -  69 

4.75 

p 

0.324 

1.606 

0.454 

*In  each  case  the  summation  is  extended  over  all  terms  for  which 


are  equal. 


APPENDIX  VII 


Values 

of  ^cR.  ,Ap, 

)(?«^^/oPj)  for 

Crystalline 

Methane 

* 

Type  of 

STRETCH 

STRETCH 

BEND 

BEND 

J 

"ij 

molecule 

z 

X,  y 

z 

X,  y 

7  - 

10 

3.13 

P 

0 

-0.175 

-0.901 

-2.114 

11 

4.06 

P 

0.666 

-0.666 

0.333 

-0.084 

20 

-  23 

4.13 

P 

-1.335 

-1.355 

-0.666 

-1.667 

32 

-  35 

4.50 

P 

0.704 

0 

0.509 

-0.764 

Ho 

-  43 

5.18 

P 

0.718 

0 

-0.297 

-1.208 

H4 

-  47 

5.54 

P 

0 

-0.963 

0.393 

-1.271 

62 

-  6? 

5.84 

P 

-2.000 

-2.000 

-2.000 

-2.000 

68 

-  69 

4.75 

P 

-0.325 

0 

-1.606 

-0.292 

1  - 

4 

2.77 

q 

-0.313 

0 

-0.178 

-0.293 

5  - 

■  6 

3.13 

q 

-0.4oJ^ 

o.4oJ^ 

-1.070 

0.267 

12 

-  19 

4.13 

q 

-2.667 

-2.667 

-3.333 

-2.333 

24 

-  27 

4.25 

q 

0.055 

0 

1.468 

-1.457 

28 

-  31 

4.50 

q 

-0.143 

0.145 

-0.914 

-1.053 

36 

-  39 

4.73 

q 

-0.429 

0 

-2.376 

-l.o4i 

48 

-  49 

5.54 

q 

0.778 

-0.778 

0.111 

-0.027 

50 

-  53 

5.60 

q 

-1.665 

0 

0.036 

-0.599 

54 

-  57 

5.73 

q 

-2.119 

0 

-0.429 

0.380 

58 

..  61 

5.79 

q 

0.970 

-0.970 

0.023 

-0.759 

In  each  case  the  trunanation  is  extended  over  all  terms  for  which  |R.  .| 
equal;  where  necessary  suius  for  p  and  q  type  molecules  are  given  separately. 


APPENDIX  VIII 


Values  of  E  for  Crystalline  Methane 

(  Units  :  cm”^  x  10°  ) 


J 

Type  of 
Molecule 

STRETCH 

X,  y,  z 

BEND 

Z' 

BEND 

X,  y 

1  -  4 

2.77 

q 

0.796 

1.114 

1.015 

5  -  10 

5.13 

p»q 

1.561 

1.175 

1.114 

11 

4.06 

p 

0.082 

0.164 

0.226 

12  -  25 

4.13 

p>q 

1.940 

1.940 

1.940 

2h  -  27 

4.25 

q 

0.818 

0.522 

0.537 

28  -  55 

4.50 

p^q 

1.218 

1.344 

1.083 

36  -  39 

4.73 

q 

0.648 

0.388 

0.962 

i+0  -  43 

5.18 

p 

0.584 

0.550 

0.444 

44  -  49 

5.54 

p^q 

0.530 

0.825 

0.815 

50  -  53 

5.60 

q 

0.352 

0.688 

0.463 

54  -  57 

5.73 

q 

0.268 

0.594 

0.457 

58  -  61 

5.79 

q 

0.4o6 

0.462 

0.501 

62  -  67 

5.84 

p 

0.684 

0.684 

0.684  ■ 

68  -  69 

4.75 

p 

0.352 

0.082 

0.326 

In  each  case  the  sumtijation  is  extended  ovei’  all  terms  for  which 
equal. 


are 


APPENDIX  rX 


Values  of  L  ‘'^R  /Sp  r^p  for  Crystalline  Methane 
J  ij  1  j 

(  Units  :  cm"^  x  10®  ) 


J 

■'n 

Type 

of 

molecule 

STRETCH 

z 

STRETCH 

X,  y 

BEND 

z 

BEND 

X,  y 

7  - 

'  10 

3.13 

P 

0 

-0.156' 

-0.350 

-0.069 

11 

4.06 

P 

-0.082 

0.082 

-0.164 

-0.144 

20 

-  23 

4.13 

P 

-0.644 

-0.61|i+ 

-0.808 

-0.566 

52 

-  55 

4.50 

P 

-0.452 

0 

-0.708 

-0.276 

4o 

-  45 

5.18 

P 

-0.396 

0 

-0.456 

-0.154 

44 

-  47 

5.54 

P 

0 

0.054 

-0.430 

-0.277 

62 

-  67 

5.84 

P 

-0.684 

-0.684 

-0.684 

-0.684 

68 

-  69 

4.75 

P 

-0.072 

0 

0.058 

-0.149 

1  - 

2.77 

a 

-0.368 

0 

-0.900 

-0.616 

5  - 

^  6 

5.15 

q. 

0.342 

-0,342 

0.150 

-0.510 

12 

-  19 

4.13 

q 

-1.292 

-1.292 

-1.128 

-1.372 

24 

-  2/ 

4.25 

q 

-0.340 

0 

-0.280 

-0.126 

28 

-  31 

4.50 

q 

0.328 

-0.328 

-0.232 

-0.562 

56 

-  39 

4.73 

q 

-0.168 

0 

-0.C60 

-0.314 

48 

o\ 

1 

5.54 

q 

-0.020 

0.020 

-0.l4o 

-0.236 

50 

-  55 

5.60 

q 

0.060 

0 

-0.484 

0.022 

54 

-  57 

5.73 

q 

0.136 

0 

-0.592 

-0.416 

58 

-  6l 

5.79 

q 

0.064 

-0.064 

-0.224 

-0.330 

In  each  case  the  summation  is  extended  over  all  tenns  for  which  IRj^I  are 
equal;  where  necessary  suras  for  p  and  q  type  molecules  arc  given  separately. 


Captlona  for  the  Figures 


Figure  1  A.  The  ahsoiTJtion  spectra  of  the  low-temperattire  phase  II  of  crystalline 
CH4  and  CD4  in  the  V3  tmd  U4  regions. 

Figure  IB.  The  absorption  spectra  in  the  V4,  region  of  2  °/o  solutions  of  both 
Isotopic  species  token  in  phase  II  of  the  solvent  species.  The  ntimbers  Inside 
the  bonds  indicate  the  line  width  in  cm”^. 


Figure  2.  Schematic  representation  of  the  crystal  structure  suggested  for 
methane  in  phase  II.  The  carbon  atoms  of  type  q  molecules  (indicated  by  a  cir¬ 
cle  at  the  carbon  atom)  are  I/2  unit  cell  above  or  below  plane  of  figure.  Two 
edges  of  each  tetrahedron  are  poreillel  to  figure  plane;  the  upper  edge  is  shown 
os  a  solid  line. 

One  face  of  the  crystallographic  unit  cell  (f.c.c.)  is  formed  by  the  five 
molecules  q,  qia,  qib.  The  primitive  unit  cell  contains  Just  the  molecules  p 
and  q. 

The  lattice  constant  in  the  f.c.c.  cell  is  h  =  A.  In  the  tetragonal 
unit  cell  the  lattice  constants  are 

c=h=5.84Aanda=b  =  h.2"^/=  =  4.13  A 

The  intermolecular  hydrogen  contacts  to  atom  pi  for  which  R  <  6a  are 
indicated  in  the  Figure.  Hydrogens  belonging  to  molecules  q  which  are  below 
plane  of  figure  are  marked  with  a  circle. 
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FIGURE  2. 


